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SUMMARY 


In the case of non-multiple eigenvalues, each of the three real eigenvalue 
extraction methods available in NASTRAN will, for a given type of normalization, 
give essentially the same eigenvectors. However, this is not so in the case of 
multiple eigenvalues. This apparent discrepancy is explained and illustrated 
by considering the example of a NASTRAN demonstration problem that has both 
multiple and non-multiple eigenvalues. 


INTRODUCTION 


The NASTRAN program provides three basic methods for real eigenvalue ex- 
traction. In each of these methods, the eigenvectors obtained can be normalized 
in three different ways (see Appendix). 

In the case of non-multiple eigenvalues, each of the three extraction me- 
thods will, for a given type of normalization, give essentially the same eigen- 
vectors. However, in the case of multiple eigenvalues, the three extraction 
methods will, in general, give different eigenvectors even though they may em- 
ploy the same type of normalization. Furthermore, in this case (of multiple 
' eigenvalues), even a given method using a given type of normalization may yield 
different eigenvectors under different conditions (e.g., different frequency 
limits on the EIGR bulk data card) [1]1. This discrepancy may seem disturbing* 
but it is explained in this paper where it is shown that the different eigen- 
vectors corresponding to multiple eigenvalues obtained by different methods and 
under different conditions have certain definite relationships among them and 


1 
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that these relationships depend on the number of arbitrary constants that are 
inherently assumed in the computation of eigenvectors. This is illustrated by 
considering the example of a NASTRAN demonstration problem that has both multi- 
ple and non-multiple eigenvalues. 


THEORETICAL BACKGROUND 


The basic eigenvalue problem solved in NASTRAN can be formulated as 



-A 



] (x> = 0 


(1) 


where K ga and M Qd are respectively the stiffness and mass matrices (both of 

which are real and symmetric) referred to the analysis set [2], A is a scalar 
quantity, and {x} is a column vector that comprises all the degrees of freedom 
in the analysis set. Non-trivial solutions for (x } in the above equation are 
possible if and only if the resultant matrix within the brackets is singular 
[3,4]. The values of A that satisfy this condition represent the required 
eigenvalues. 


Let A^ be one of the eigenvalues of Equation (1). The eigenvector cor- 
responding to this eigenvalue can be obtained from the equation by substituting 
* = A^ and solving for (x> . We thus have 

l*aa M aal { *><t = 0 < 2 > 

where {x}^ is the required eigenvector. 

Let n be the order of the problem and let x^, x 2 , x^, x p be the com 

ponents of (x). Equation (2) can then be rewritten as 

n 

E C ,, x. = 0 , i = 1,2,3, n (3) 

j=l J J 


where C.^ are constants that depend on K aa , M aa , and A^. 

Equation (3) represents a system of n linear equations in n unknowns. 
However, not all of these n equations are independent. The exact number of in- 
dependent equations depends on the multiplicity of the eigenvalue H- ^ s 
(1<. sin) is the multiplicity of A$, , then it can be shown that the rank of the 
resultant matrix in Equation (2) is (n-s) [3,4]. The number of independent 
equations in (3) is also, therefore, (n-s). Their solution thus involves s 
arbitrary constants. The total number of eigenvectors available is therefore 

equal to oo s . However, it can be shown that the number of linearly independent 
eigenvectors is only s [4]. 
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Depending on the multiplicity s, it is useful to distinguish between the 
■oil owing two cases: 

:ase (a) Non-multiple Eigenvalues (s = 1) 

In this case, (n - 1) of the n equations in (3) are independent and their 
;olution involves a single arbitrary constant. An infinite number of eigen- 
/ectors are thus available. However, the important thing to note here is that 
:he relative values of the eigenvector components remain invariant . 

)ase (b) Multiple Eigenvalues (1 <s<ji) 

In this case, the number of independent equations in (3) is equal to (n - s) 
md their solution involves s arbitrary constants. In contrast to the case of 
^on-multiple eigenvalues, the most important thing to note in this case is that 
the relative values of the eigenvector components are not invariant, but depend 
j>n the relative values of the involved constants themselves. Also, it is clear 
pat the solution space in this case is much larger than in the case of non- 
pul tiple eigenvalues. 


EIGENVECTOR COMPUTATION IN NASTRAN 


The exact procedure for eigenvector computation in NASTRAN depends on the 
method used for eigenvalue extraction and is described in detail in the Theo- 
retical Manual [2]. Thus, in the Inverse power method, trial eigenvectors are 
Assumed and iterated until convergence occurs. In the Determinant method, 
eigenvectors are computed by the method of backward substitution after each of 
the corresponding eigenvalues has been calculated. In the Givens method, all 
the eigenvalues are first obtained and then the required number of eigenvectors 
is computed by repeated use of backward substitution. 

The orthogonality of eigenvectors of closely spaced eigenvalues is guar- 
anteed by the procedure employed in the Inverse power method. NASTRAN employs 
the Gram-Schmidt orthogonal izati on procedure [4] for the purpose in the case of 
ithe Determinant and Givens methods. 

Whatever the method employed, the eigenvectors obtained all exhibit the 
(characteristics described in the previous section and also involve the inherent 
use of one or more arbitrary constants. In general, for a given mode, the in- 
volved constants for the computation of the corresponding eigenvector will not 
be the same in all the three methods. It should also be noted that the user has 
|no control whatever on the selection or the choice of these constants. 
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EXAMPLE 


To illustrate the points made in the paper, we consider the problem of the 
vibration of a square plate with hinged supports at all the four edges. (This 
is the same as NASTRAN Demonstration Problem 3-1. See Reference 5 for details.] 

The finite element model employed is shown in Figure 1. Only half the 
plate has been modelled and symmetric boundary conditions used on the center-1 ir 
in order to reduce the order of the problem as well as the bandwidth. 

The parameters of the model are as follows: 

Length = l = 508.0 mm (20.0 in.) 

Width = w = 508.0 mm (20.0 in.) 

Thickness = t = 25.4 mm (1.0 in.) 

Young's modulus = E = 2.06843 x lO^ 1 N/m 2 (3.0 x 10^ psi) 

Poisson's ratio = v = 0.3 1 

Mass density = p = 2.20197 x 10^ Kg/m^ (206.04393 lb.-sec 2 /in.^)i 

The eigenvalues and eigenvectors for the first six modes of the model were , 
obtained by the Inverse power and Determinant methods. 2 For each method, four 
separate runs were made with different eigenvalue extraction data. The detailed 
data for all the eight runs are given in Table 1. Since the number of degrees 
of freedom in the analysis set (or a-set) in this problem is as large as 590 
and only the first few modes were required for the purpose of the present paper, 
the Givens method was considered unsuitable. 

The eigenvalues and natural frequencies obtained in the eight runs (as well 
as the corresponding theoretical values [6]) are presented in Table 2. It can 
be seen from this that the third and fourth modes together represent an eigen- 
value of multiplicity two, that is, a double root. The other modes yield non- 
multiple eigenvalues. 3 

Tables 3 through 7 present, for the first six modes, the eigenvector com- 
ponents corresponding to the vertical displacements at points 111 through 121 
(along the line AB in Figure 1). The results for the third and fourth modes, 
which represent a multiple eigenvalue, are included together in Table 5. 


^ The runs were made on the CDC computer using an improved version of the 15.5 
level of NASTRAN. 

It is true that no conclusions can be drawn from Table 2 regarding the multi- 
plicity of the eigenvalue for the sixth mode because, in general, at least 
(m + 1) eigenvalues must be obtained before one can examine the multiplicity 
of the mth eigenvalue. However, in the present problem, the non-multiplicity 
of the sixth eigenvalue was confirmed from additional runs with more than six 
eigenvalues. 
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An examination of Table 3 (corresponding to the non-multiple eigenvalue of 
the first mode) reveals that the eigenvector components for all the runs are 
(essentially the same. (Actually, the results for all the runs except Run 3 are 
identical. The extremely small differences between the results of Run 3 and 
.those of others are so inconsequential as to be practically meaningless.) The 
results thus show that they are essentially unaffected by differences in the 
eigenvalue extraction data. 

i Tables 4, 6 and 7 (corresponding to the non-multiple eigenvalues of the 
second, fifth and sixth modes respectively) need some explanation. Note that 
all the eigenvector components given in these three tables are extremely small 
in magnitude. Actually, the corresponding theoretical values are all zero [6]. 
For all practical purposes, the differences among the various runs in these 
tables can, therefore, be considered as quite insignificant. 

Table 5 (corresponding to the multiple eigenvalue of the third and fourth 
modes) is interesting. It can be seen that the results given by the runs em- 
ploying the Determinant method are essentially the same and that they are 
[essentially unaffected by differences in the eigenvalue extraction data. The 
same thing is not wholly true for the runs employing the Inverse power method 
jbecause slight variations are noticeable among some of them. (In this con- 
mection, it should be emphasized that the association of a particular eigen- 
vector with a particular mode number in this table is irrelevant because the 
numbering of the modes for this multiple eigenvalue is completely arbitrary.) 
The most important point to note in this table, however, is that the eigen- 
vectors given by the Inverse power method are completely different from those 
given by the Determinant method. At first sight, this discrepancy in the 
results may seem disturbing. However, it can be explained from a theoretical 
viewpoint and the discrepancy may be shown to be apparent and not real. 

Now, since the third and fourth modes involve a multiple eigenvalue of 
order two, it follows from the theory presented earlier that their eigenvectors 
involve two arbitrary constants in them and that the relative values of the 
eigenvector components depend upon the relative values of these two constants. 
The differences among the various eigenvectors shown in Table 5 are therefore 
clearly due to the fact that they represent different relative values of the 
two involved constants. 

The theory also shows that, in the case of a multiple eigenvalue of order 
; two, the total number of eigenvectors possible is doubly infinite. The solu- 
tion space is thus much larger than in the case of non-multiple eigenvalues. 

The theory further shows that the number of independent eigenvectors in this 
case is only two. This means that every eigenvector in the solution field can 
be obtained by the linear combination of any two eigenvectors that are them- 
selves linearly independent. This can be shown to be true in the case of the 
eigenvectors given in Table 5. 
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Let {x} i • represent any eigenvector in Table 5 where i denotes the mode 

number and j denotes the run number. It can be seen from the table that there 
are several eigenvectors that are identical and, therefore, linearly dependent. 
One such set, for instance, is given by the eigenvectors (x ) 3 -j and (x > 42 * 

Another such set is given by {x> 32 and (x> 43 . Yet another such example is the 

set comprising {x} 35 , {x) 36 , {x> 37 and (x> 38 . A fourth such set is given by 

{x}^, {x> 47 and {x) 4 g. 


As an example, consider the eigenvectors {x> 3g and {x> 4g in Table 5 which | 

are clearly linearly independent. It can be shown that every eigenvector in 
Table 5 can be obtained as a suitable linear combination of these two eigen- 
vectors. Simple algebra shows that, for instance, the eigenvector l*} 32 in j 

Table 5 can be obtained from these two independent eigenvectors by the following ( 
(approximate) relationship: 


(x) 32 = 0.9850782(x) 36 - 0 . 1 91 71 1 6{ x> 46 (4) j 

Similarly, it can be shown that the following (approximate) relationship 
can be used to obtain the eigenvector (x) 34 from (x} 36 and {x> 4g : 1 

(x) 34 = 0.9850979{x> 36 - 0.1914862{x) 4g (5) 


In a like manner, similar linear relationships can be shown to exist be- 
tween any eigenvector and any two linearly independent eigenvectors in Table 5. 

The above discussion and relationships clearly show that the discrepancy 
in the results presented in Table 5 is thus only apparent and not real. 

It is interesting to note that the plate problem considered in this 
example yields multiple roots at many higher modes also. Thus, for instance, 
modes 9 and 10 represent a double root. So do modes 12 and 13. Results simi- 
lar to the above can be expected to occur in these cases also. 


CONCLUSIONS 


NASTRAN provides three basic methods for real eigenvalue extraction. In 
the case of non-multiple eigenvalues, each of these three methods will, for a 
given type of normalization, give essentially the same eigenvectors. However, 
this is not so in the case of multiple eigenvalues. This apparent discrepancy 
has been explained and illustrated by considering the example of a NASTRAN de- 
monstration problem that has both multiple and non-multiple eigenvalues. 
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APPENDIX 


The NASTRAN program provides three basic methods for real eigenvalue ex- 
.raction. These are: 

) The Inverse power method, 

) The Determinant method, and, 

j) The Givens (or Tridiagonal) method. 

The Inverse power and Determinant methods are "tracking" methods and as 
efficient in those cases where only a few of many eigenvalues are de- 
ired; On the other hand, the Givens method is a transformation method that is 
(fncient only in those cases where all eigenvalues or a high proportion of all 
jigenvalues are required [2], 

| In each of the three methods, the eigenvectors obtained can be normalized 
n any one of three different ways. The three types of normalization (N0RMS) 
vail able are [1]: 

) normalizing to unit value of the generalized mass (MASS), 

I 

' to unit value of the largest component in the analysis set 

(MAX), and, 

) normalizing to unit value of a specified component in the analysis set 
(P0INT). 


i 
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Table 1. Eigenvalue Extraction Data* for the NASTRAN Runs 
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See Reference 1 for details 


















Table 2. Eigenvalues and Natural Frequencies 


Mode 

no. 

NASTRAN results 

(Runs 1 through 8 — See Table 1) 

Theoretical values , 

(Reference 6) 

Eigenvalues 

(rad/sec) 2 

Natural frequencies 
(Hz) 

Eigenvalues 

(rad/sec) 2 

Natural frequencies 
(Hz) 

l a 

3.237408E+01 

9.055634E-01 

3.246970E+01 

| 

9.068997E-01 1 

— 1 

2 a 

2.022407E+02 

2.263364E+00 

2.029356E+02 

2.267249E+00 

3 a 

8.111 597E+02 

4.532870E+00 

8.117425E+02 

4.534499E+00 

4 b 

8.111 597E+02 

4.532870E+00 

8.11 7425E+02 

4.534499E+00 

5 C 

1.352052E+03 

5.852169E+00 

1.371845E+03 

5.894848E+00 

6 d 

2.355330E+03 

7.724066E+00 

2.345936E+03 

7.708648E+00 


a Obtained by all runs, 
k Obtained by all runs except Runs 1 and 5. 
c Obtained by Runs 3, 4, 7 and 8. 

^ Obtained only by Runs 4 and 8. 
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Eigenvalue = 3.237408E+01 (rad/secr Natural frequency = 9.055634E-01 Hz 
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The eigenvalue analyses were done by using the British system of units. However, the eigenvector components obtained represent 
relative values and thus are independent of the system of units employed. 




Eigenvalue = 8.111597E+Q2 {rad/ sec) Natural frequency = 4.532870E+00 Hz 
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See note given under Table 



Table 6. Eigenvector Components* Corresponding to Vertical Displacements for Modj 
Eigenvalue = 1 .352052E+03 (rad/sec) 2 Natural frequency = 5.852169E+00 Hz 


Grid point 
(see Figure 1) 

Run 3 

Run 4 

Run 7 

1 

Run 8 

111 

-5 . 69991 2E-07 

-1 .843495E-06 

-2.1 83520E-1 2 

—7 . 436551 E- 1 2 


112 

-6.080192E-07 

-1 . 758002E-06 

-2.089282E-12 

-6 . 985661 E- 1 2 

113 

-7.11 3474E-07 

-1 . 51 731 8E-06 

-1 .771 980E-1 2 

| 

-5.644270E-1 2 

1 

1 

114 

-8.5011 40E-07 

-1 .165451 E-06 

-1 . 1 35621 E-l 2 

—3 . 503232 E- 1 2 

1 

115 

-9.82041 1 E-07 

-7.652404E-07 

-2.059582E-13 

-8.178176E-13 


116 

-1 .061947E-06 

-3.849968E-07 

7.1 25514E-13 

1 . 7461 65E- 1 2 

117 

-1 .052001 E-06 

-8.408051 E-08 

1.267204E-12 

3. 572532E-1 2 , 

118 

-9 . 304697E-07 

9.946076E-08 

1.520344E-12 

4.248741 E-l 2 

119 

-6.970452E-07 

1 . 565069E-07 

1.396781 E-l 2 

! 

3.669778E-1 2 

120 

-3.735721E-07 

1 .080931 E-07 

8.6581 39E- 1 3 

2.11121 2E- 1 2 

121 

0.0 

0.0 

0.0 

0.0 


Inverse 

power method 

Determinant method 


See note given under Table 3. 
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Table 7. Eigenvector Components* Corresponding to Vertical Displacements for Mode 6 
Eigenvalue = 2.355330 (rad/sec)^ Natural frequency = 7.724066E+00 Hz 


Grid point 
(see Figure 1) 

Run 4 

Run 8 

111 

-4.899823E-10 

-6.518056E-1 2 

112 

-4.359024E-10 

-5.832685E-1 2 

113 

-2.854791E-10 

-3.934395E-12 

114 

-7.157563E-11 

-1 - 2521 90E-1 2 

115 

1 . 590946E-1 0 

1 .593764E-1 2 

116 

3. 561466E-10 

3.982927E-1 2 

117 

4.764730E-10 

5.41 4795E- 1 2 

118 

4.937036E-10 

5.61 2997E- 1 2 

119 

4.039232E-1 0 

4. 588831 E- 12 

120 

2.265359E-10 

2 . 571 805E— 1 2 

121 

0.0 

0.0 


Inverse power method 

Determinant method 


See note given under Table 3. 
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